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Abstract: Congenital anomalies are the main causes of preterm and neonatal mortality and 
morbidity. We investigated the association between congenital anomalies and mothers' 
exposure to air pollution during pregnancy by combining risk estimates for a variety of air 
pollutants (SO2, NO2, PM10, PM2.5, CO and O3) and anomaly defect outcomes. 
Seventeen articles were included in the systematic review and thirteen studies were taken 
into account in the meta-analysis. Combined estimated were calculated separately 
according to whether the exposure metric was continuous or categorical. Only one 
significant combination was; NO2 concentrations were significantly associated with 
coarctation of the aorta (OR = 1.20 per 10 ppb, 95% CI, (1.02, 1.41)). This finding could 
stem from strong heterogeneity in study designs. Improved exposure assessment methods, 
in particular more accurate spatial measurements or modeling, standardized definition of 
cases and of better control of confounders are highly recommended for future congenital 
anomalies research in this area. 
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1. Introduction 

Today, air pollution is a major public health concern. Due to industrial emissions, urbanization and 
transport of goods and people by fuel-engine motor vehicles, air pollution affects everyone in developed 
and developing countries. Common pollutants such as nitrogen oxide (N0 2 ), sulfur dioxide (S0 2 ), 
particulate matter (PM) and carbon monoxide (CO) have been shown to be associated with several adverse 
health events such asthma attacks and incidence [1,2], chronic obstructive pulmonary diseases [3], 
cardio and cerebrovascular conditions and lung cancer [4-6], and have been reported to reduce human 
fertility [7]. Environmental nuisances affect all age groups. However, fetus development, newborns and 
infants are recognized to be more particularly vulnerable to air pollution [8-10]. Exposure of pregnant 
women is linked with fetal growth retardation, low birth weight [11], preterm birth and neonatal 
mortality [12]. Maternal exposure to air pollution may be also related to congenital anomalies. 
However, the evidence is still weak due to the paucity of epidemiological studies. 

Congenital anomalies are recognized to be a major risk factor of stillbirth and of neonatal and infant 
mortality. Worldwide, an estimated 10% deaths under five-years-old children are caused by congenital 
anomalies [13]. European Surveillance of Congenital Anomalies (EUROCAT), an organization of 
population-based registries for the surveillance of congenital anomalies in Europe, recorded the 
perinatal death rate 9.3 per 10,000 births of all congenital anomaly between 2008 and 2012. Among 
them, 23.7% was due to chromosomal anomalies, 22.6% to congenital heart defects, and 17.2% to 
nervous system anomalies [14]. About half of all major congenital malformations are of unclear 
etiology and are suggested to have multifactorial causes, including environmental exposures [15]. 

In the past decade, the number of studies investigating the association between congenital 
anomalies and air pollution has increased. The potential impact of environmental exposures to 
congenital anomalies has been recently reviewed [16,17]. The studies concluded that exposure to NO2, 
SO2 and PM2.5 increased significantly the risk of congenital heart diseases. Since the last 2011 review, 
several new studies have been published [18-23]. In this setting, updating the literature synthesis may 
improve our understanding of the relationship between air pollution and congenital anomalies and also 
of the biological process through which air pollution could lead, directly or indirectly, to these 
outcomes. We, therefore, conducted a meta-analysis to assess the association between air pollution and 
the risk of congenital anomalies. We explored whether some subtype of anomalies could be particularly 
concerned by air pollution. Finally, we discussed hypotheses explaining the different routes by which 
air pollution might increase the risk of congenital anomalies. 

2. Methods 

The methodology adopted has been described in detail in the previous published review [16]. 
Meta-analyses were conducted for a minimum number of four individual studies. We separated the data 
set into two categories according to the metric used for exposure assessment (continuous or categorical). 
Combined odds ratios were computed in order to contrast the highest with the lowest quartiles of exposure 
when the individual studies reported categorical exposure metrics. When exposure was expressed as 
a continuous variable, if quantitative descriptors for air pollutants were available, we also converted 
effect estimates into ORs contrasting the highest versus the lowest exposure categories [24]. 
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Otherwise, risk estimates which had been calculated from continuous exposure metrics were expressed 
as unit odds ratios, corresponding to an increase of 1 microgram per cubic meter (ug/m ) for sulfur 
dioxide (SO2), 10 ug/m 3 for particulate matter with diameter less than 10 urn (PM10) and to an increase 
of 1 part per million (ppm) for carbon oxide (CO), and 5 part per billion (ppb) for ozone (O3) and 10 ppb 
nitrogen dioxide (NO2). Key features and definitions of exposure of each study are detailed in Table 1. 

2.1. Search Methods 

A literature search was conducted in the PUBMED database in order to select articles published 
between January 2011 and January 2014. The search strategy followed the PRISMA guidelines [25]. 
The keywords used for this review were (air pollution OR traffic pollution OR outdoor pollution) 
AND congenital anomalies. We also used the terms "traffic pollution" and "outdoor air pollution". 
Searches were restricted to English-language articles. No restriction was put on the geographical 
location. Abstracts of all studies were then screened manually and excluded if they were not performed 
on human populations and did not present original data (review articles). Full manuscripts were 
checked thoroughly. Seven studies were published after 2011; one was not included because its main 
issue did not deal with the association between congenital anomalies and air pollutants [26]. 
We also included the eleven articles used in the previous literature synthesis published in 2011. 
Overall, seventeen articles were included. 

2.2. Data Extraction 

We selected measures from the adjusted models presented in each study. Odds ratios and similar 
metrics relating outcomes and pollutants were extracted. For cohort studies, we used risk ratios since 
the two ratios give equivalent results when the outcome is rare. In addition, the period of exposure 
during the pregnancy has been taken into account since the pregnancy weeks 3-8 constitute the critical 
window of exposure for embryogenesis; later exposures may not contribute to the etiology of major 
congenital anomalies [22]. 

2.3. Meta-Analysis 

Heterogeneity was assessed for pollutant-congenital anomaly outcomes by using the Cochrane 
Q-test. Fixed effect models were used when the result of the Q-test gave a heterogeneity p-value 
higher than 0.1. In contrast, random effects models were used for p-values lower than 0.1. 
Following Higgins et al. [27], a low heterogeneity was determined for I between 25%-50%, 
moderate between 50%-75%, and high for >75%, where I 2 is defined as the percentage of variation 
attributable to heterogeneity. Forest plots were generated to illustrate the combined risk estimates. 
Statistical analysis was performed using STATA 1 1 (TX, USA). 
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Table 1. Overview of studies included in the systematic review. 



Study 



Location Period Study Design Congenital Categories Exposure Assessment Exposure Variable Air Pollutants 



Results 



Confounders 



Gianicolo 
et al. 2014 
[23] 



Case-control, 
Brindisi, Italy 2000-2010 individual 

matching 



Daily average 
concentration of 
Congenital heart defects, pollutants measured by 3 



atrial septal defects 



monitoring stations and 
performed for week 3-8 
of gestation 



Continuous and 
categorical 



SO, and TSP 



Schembari 
etal. 2013 
[22] 



Barcelona, 
Spain 



1994-2006 



Case-control, 
no matching 



Congenital heart defects, 
neural tube defects, 
respiratory system 
defects, orofacial clefts, 
digestive system defects, 
abdominal wall 



Daily spatio-temperal 
exposure estimates over 
week 3-8 of pregnancy 



Continuous 



N0 2 , NO„ PM 1( 
PM, PIL™ 



Exposure to 90th 
percentile of SO2 
increased risk of CHD 
(p = 0.01) and VSD 
(p < 0.05) 

Statistically significant 
associations (p = 0.05) 

between N0 2 and 
coarctation of the aorta 
and digestive system 
defects, and between 
)arse and abdon 
wall defects 



No adjusted 
confounders; cases 
and controls were 
matching for gender, 

socio-economic 
deprivation and the 
year of pregnancy 



Maternal age, 
socio-economic 
status, year of birth, 
conception season 



Padula etal, 
2013 
[20] 



California, 
USA 



1997-2006 



Case-control, 
no matching 



Anotia/microtia, 
anorectal atresia/stenosis, 
craniosynostosis, hypospadias 
degree, diaphragmatic hernia, 
transverse limb deficiency, 
intestinal atresia/stenosis, 
amniotic band syndrome, 
limb body wall complex, 
hydrocephaly, 
longitudinal limb deficiency, 
esophageal atresia 



Residence-based 
assignments around 
stations, with daily 
average values during 
first two months of were 
collected; a maximum 
interpolation radius 
of 50 km was used 



Categorical 



N0 2 , NO, CO, 
O3, PM 10 , PM2.5 



No significant 
association had been 



revealed 



Maternal ethnicity, 
education, and early 
prenatal vitamin use 
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Table 1. Cont. 



Study 


Location 


Period Study Design 


Congenital Categories 


Exposure Assessment 


Exposure 
Variable 


Air Pollutants 


Results 


Confounders 










Residence-based 


















assignments around stations, 










Padula 
etal. 2013 
[19] 


California, 
USA 


Case-control, 

1997-2006 

no matching 


Congenital heart diseases 
groups (27 subtypes) 


with daily average values 
during first two months of 
were collected; a maximum 
interpolation radius 
of 50 km was used 


Categorical 


NO, NO z , PMm, 
PM 2 .5, CO, O a 


No significant 
association had been 
revealed 


Maternal 
ethnicity, education, 
and vitamin use 


















Infant's sex, 


















plurality, religion, 








Multiple congenital heart, 


Weekly means of exposures 








maternal age, 


Agay-Shay 

K et al. 
2013 [21] 


Tel-Aviv, 
Israel 


Case-control, 

2000-2006 

no matching 


atrial and atrial 
septal defects, 
isolated ventricular 
septal defects 


during pregnancy week 3—8 
according to the distance from 
stations to each maternal 
address 


Continuous 


N0 2 , S0 2 , 
PMio, PM 2 .5, 
CO, 0 3 


No significant 
association had been 
revealed 


maternal and paternal 
marital status, 
maternal and 
paternal origin, 

natprnfll aop and tbp 

season of conception 










Residence-based assignments 










Padula 
et al. 2013 
[18] 


California, 
USA 


Case-control, 

1997-2006 

no matching 


Neural tubes defects (spina 

bifida and anencephaly), 
orofacial clefts, gastroschisis 


around stations, with daily 
average values during first two 

months of were collected; 
a maximum interpolation radius 


Categorical 


N0 2 , NO, 
PM10, PM 2 5 , 
CO, o 3 


No significant 
association had been 
revealed 


Maternal ethnicity, 
education, and 
vitamin use 



of 50 km was used 
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Table 1. Cont. 



Study 



Location 



Period Study Design Congenital Categories Exposure Assessment Exposure Variable Air Pollutants 



Results 



Confounders 











An association between 








Cardiac chambers and 




NO2 and congenital heart 








connection, cardiac septa, 




diseases, ventricular septal 


Socio-economic 












pulmonary and tricupid 




defect, cardiac septa 


status, degree of 












valves, aortic and mitral 


Weekly average of 


malformations and 


Dadvand 


Case-control, 




urbanity, and season 




Northeast of 


valves, great arteries and 


pollutants at nearest 


S0 2 , N0 2 , CO, tetralogy of Fallot; and 


et al. 2011 


1993-2003 frequency 


Continuous 


of conception; 




UK 


veins, atrial septal defect, 


monitors to maternal 


PM10, 0 3 CO exposure to 


[28] 


matching 




cases and controls 


coarctation of aorta, 


residential location 


ventricular septal defect, 












were matching for 






pulmonary valve stenosis, 




cardiac septa 








the year of birth 






tetralogy of Fallot, 




malformations and with 






ventricular septal defect 




congenital pulmonary 












valve stenosis 





Dadvand 
etal. 2011 
[29] 



Northeast of 
UK 



1985-1996 



Case-control, 
frequency 
matching 



Coarctation of aorta, 
tetralogy of fallot, congenital 
pulmonary value stenosis, 

atrial septal defect, 
ventricular septal defect, 
congenital cardiac chambers 
and connections, congenital 
cardia septa, congenital 
pulmonary and tricuspid 
valves, congenital aortic and 
mitral valves, congenital great 
arteries and veins 



Weekly exposure 
levels by two stage 
spaiotemporal modeling 
at each maternal place 
of residence 



Categorical and 
continuous 



so 2 , 

black smoke 



Birth year, 

An association between socio-economic 
maternal exposure to status, infant sex, 
black smoke and cardiac season of conception, 

chambers and and degree of 

connections (only when urbanity; cases and 
using exposure as a controls were 

continuous variable) matching for the 
year of birth 
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Table 1. Cont. 



Congenital Exposure 
Study Location Period Study Design Exposure Assessment Air Pollutants Results Confounders 

Categories Variable 

















Maternal age, 
















ethnicity, smoking 


Marshall 
et al. 2010 
[30] 


Case-control, 

New Jersey, 

1998-2003 frequency 

US 

matching 


Cleft lip with cleft palate, 
cleft palate 


Average concentration 
of exposures at nearest 
monitor stations 
(13-20 km) 


Categorical 


S0 2 , N0 2 , CO, 
PMio, PM 2 .5, 0 3 


No significant 
association had been 
revealed 


and alcohol, and 
season of conception; 
cases and controls 
were matching with 
maternal residence at 
birth 






Anomalies of cardiac 
















chambers, transposition of 












Dolk 
etal. 2010 
[31] 


Wessex, 
North West 
Thamas, 

1991-1999 Cohort 

Oxford and 
Northern of 
UK 


great vessels, 
malformations of cardiac 

septa, atrioventricular 
septal defects, tetralogy of 
Fallot, malformations of 
valves, hypoplastice left 


Annual mean 
exposure at census 
level in 1996 


Continuous 


S0 2 , N0 2 , 
PM 10 


A significant 
association between 
S0 2 and tetralogy of 
Fallot, and between 
PM 10 and omphalocele 


Maternal age, 
socio-economic 
deprivation 



heart syndrome, great 



arteries and veins, 
coarctation of aorta 
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Table 1. Cont. 



Study Location Period Study Design Congenital Categories Exposure Assessment Exposure Variable Air Pollutants 



Results 



Confounders 



Hansen 
et ai. 2009 
[32] 



Brisbane, 
Australia 



1998-2004 



Case-control, 
individual 
matching 



Aortic artery and valve 

defects, atrial septal 
defects, pulmonary artery 

and valve defects, 
ventricular septal defects, 

conotruncal defects, 
endocardial cushion and 
mitral valve defects, cleft 
lip, cleft palate, cleft lip 
with cleft palate 



Daily average exposures 
at 1 8 monitors with the 
majority located within 
a 30 km radius of city 



Continuous 



SO2, NO, CO, No significant association 
PM 10 and O3 had been revealed 



Infant sex, birth order, 
season of birth, 
maternal age, 
education, alcohol, 
and body mass index; 
cases and controls 
were matching with 
mother's age, 
marital status, 
number of previous 
pregnancies, 
month of LMP, 
area-level SES, 
and distance to 
pollution monitor 



Rankin 
et at, 2009 
[33] 



Nervous system, 
congenital heart defects, 
atrio ventricular septal 
defects, tetralogy of fallot, 

hypoplastic left heart, 
coarctation of aorta, patent Daily average exposures 



Northern 
region, UK 



1985-1990 



Case-control, 
no matching 



ductus arteriosus, 
ventricular septal defect, 
respiratory tract, cleft lip 
and palate, eye, ear, face 

and neck, digestive 
system, internal urogenital 
system, musculoskeletal, 
miscellaneous 



during the first trimester 
from monitors within 10 
km of maternal residence 



Continuous and 
categorical 



S02, black smoke 



A significant association 
between black smoke and 
nervous system 
anomalies. 



Birth weight, infant 
sex, and maternal 
deprivation 
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Table 1. Cont. 



Study Location Period Study Design Congenital Categories 



Exposure Assessment 



Exposure 
Variable 



Air Pollutants 



Results 



Confounders 



Strickland 
et al. 2009 
[34] 



Atlanta, 
USA 



1986-2003 



Cohort 



Tansposition of the great 
arteries, tetralogy of fallot, 
pulmonary stenosis and valvar, 
patent ductus arteriosus, 
hypoplastic left heart , 
coarctation of aorta, 
atrial septal defect, secundum, 
ventricular septal defect, 
muscular, ventricular septal 
defect, permimemnranous, 
conotruncal defect, 
Left ventricular outflow tract 
defect, right ventricular outflow 
tract defect 



Average of daily 
concentration from 
one central monitoring 
station 



A significant association 
S0 2 , N0 2 , CO, No adjusted 

Continuous between PM 10 and patent 

PMio, and O3 variables 
ductus arteriosus 



Hwang 
et al. 2008 
[24] 



Taiwan 2001-2003 



Case-control, 
no matching 



Cleft lip 



Monthly average of 
exposures at 72 stations 
by using inverse 
distance weighting 
method during the 
first trimester 



Continuous 



S0 2 , NO z , CO, 
PM 10 and 0 3 



A significant association 
for first and second month 
O3 exposure 



Maternal age, infant 
sex, plurality § , 
gestational age, 
population density, 
and season of 
conception 
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Table 1. Cont. 



Study 



Location 



Period 



Study 
Design 



Congenital 
Categories 



Exposure Assessment 



Exposure 
Variable 



Air Pollutants 



Results 



Confounders 



















Infant sex, birth order, 


Kim 
et al. 2007 
[35] 


Seoul, 
Korea 


2001-2004 Birth cohort 


Not specified 


Residence based average 
exposure levels at each 
trimester from nearest 
monitoring stations 


Continuous and 
categorical 


PM 10 


Congenital anomalies 
were influenced by 
exposure to PM 10 


season of birth, maternal 
age, maternal and paternal 
education, alcohol, body 
mass index and maternal 
weight before delivery 



Giloba 
et al. 2005 
[36] 



Texas, 
USA 



1997-2000 



Case-control, 
frequency 
matching 



Aortic artery and 
valve defects, atrial 
septal defects, 
pulmonary artery and 
valve defects, 
ventricular septal 
defects, conotruncal 
defects, endo-cardial 

cushin and mitral 
valve defects, cleft lip 
with cleft palate, cleft 
palate 



Average of daily 
measurements based on 
the first closest monitor 

(median distance 
8.6-14.2 km) 



Categorical 



S0 2 , N0 2 , CO, 
PM 10 and Q 3 



A significant association 
between exposure to S0 2 
and VSD(p< 0.0001), 
CO and tetralogy of fallot 



Maternal age, ethnicity, 
education, marital status, 

illness, tobacco use, 
season § of conception, 
plurality, parity, infant sex, 
prenatal care, 



(p < 0.0017), PM 10 and ASD and gravidity; cases and 
(p < 0.0001), S0 2 with ASD control were matching with 
(0.0017) vital status, year, 

maternal county of 
residence at delivery 
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Table 1. Cont. 



Study Location Period Study Design Congenital Categories Exposure Assessment Exposure Variable Air Pollutants Results Confounders 









Aortic defects, 










Maternal age, ethnicity, 








Pulmonary valve, 










education, marital status, 








Conotruncal defects, 










illness, tobacco use, 


Ritz 
et al. 2002 
[9] 


California, 
USA 


Case-control, 

1987-1993 

no matching 


Ventricular septal defects, 
Multiple cardiac or 
cleft defect, Syndrome 
with cardiac or cleft 
defect, Isolated cleft 
palate, Isolated cleft lip 
with/without palate 


24 h average 
measurements every 
6 days over duration of 
pregnancy 


Continuous and 
categorical 


CO, o 3 


No significant 
association had been 
revealed 


season of conception, 

plurality § , parity, 
decade of infant's birth, 
infant sex, access to 
prenatal care, 
time since last pregnancy 
and birth type 



Notes: S0 2 , sodium dioxide; NO, nitrogen oxide; CO, carbon oxide; PM 10 , particulate meter with diameter < 10 um; PM 25 , particulate meter with diameter < 2.5 urn; 0 3 , ozone; 
TSP, total suspended particulate; CHD, congenital heart defects; ASD, atrial septal defects. s Plurality stands for multiple births (including stillbirth) after 
one pregnancy. 
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3. Results 

Seventeen articles met the inclusion criteria for the systematic review, and the characteristics of 
each study are shown in Table 1. Three articles published in 2013 were conducted by the same 
research group [18-20]. We considered them as three independent studies because each dealt with 
different outcomes. Another team also split results into two articles [28,29]. Seven studies were 
conducted in the United States [9,18-20,30,34,36], four in the United Kingdom [28,29,31,33] and only 
one in Australia [32], in Israel [21], in Italy [23], in South Korea [35], in Spain [22] and in Taiwan [24]. 
Overall, twelve studies were case-control studies and five were cohort studies. One did not specify 
congenital subgroups and diagnosed only birth defects (n = 14); it was not included in the 
meta analysis [35]. 

Most studies used a population-based case-control design, selecting cases from clinical or autopsy 
reports. Controls were randomly selected from birth registries. Cases included live birth, stillbirth or 
termination of pregnancy after a congenital anomaly diagnosis. One study collected only newborns with 
diagnosis of congenital anomalies [23]. Two studies focused specifically on orofacial defects [24,30] 
and therefore included only cases with a diagnosis of cleft lip with or without cleft palate. Definitions of 
cases with single or multiple congenital defects and criteria for splitting them into sub-groups varied 
across studies. Cases definition used in the studies published after 2010 were mainly coded according 
to the International Classification of Diseases, version 9 or 10 (ICD codes). Cardiovascular anomalies 
were the most frequently investigated defects, followed by neuron system. 

Mother residence concentration estimates from air quality monitoring networks were frequently 
used for exposure assessment. Average concentrations of pollutants during the first or the first two 
months of pregnancy were calculated from the nearest monitors. The distance from monitors to 
maternal residence varied among studies, from 10 km to a maximum of 50 km. Classical air pollutants 
were measured, i.e., SO2, NO2, PM10, PM2.5, CO and O3 in most studies. Total suspended particulate 
(TSP), NO, NO x and black smoke were assessed in a small number of studies. Pollutant concentration 
distributions are shown in Table 2. 

We conducted the meta-analysis for 21 combinations of air pollutants and congenital anomalies 
when at least four studies were available for the same combination. Heterogeneity tests (the Q-test) 
indicated four combinations with high values of I , for which random effects models were applied. 
Heterogeneity varied between 0% and 93.4%, which indicated that the measurement methods, sample 
property and characteristics varied among and within different congenital groups. 
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Table 2. Exposure distribution in studies included in the meta analysis. 



Studies 


Case 


Control 


Total 


Congenital Defects Air Pollutants 


Congenital Defects Air Pollutants 


Gianicolo et al. 


CHD S0 2 , mean 2.9 (ig/m 


CHD S0 2 , mean 2.8 (ig/m 




2014 [23] 


VSD S0 2 , mean 3.2 |ig/m 3 


VSD SO z , mean 2.8 ng/m 3 




Schembari et al. 


IQR: N0 2 , 12.7 ug/m 3 


N0 2 , IQR 11.8 (ig/m 3 




2013 [22] 


PM 10 , 2.8 |ig/m 3 


PMio, IQR 3.0 ng/m 3 










Minimum, median, maxmum: 








S0 2 0.33 ppb, 2.1 ppb, 51.4 ppb; 


Agay-Shay et al. 






CO, 0.15 ppm, 0.9 ppm, 13.5 ppm; 


2013 [21] 






N0 2 0.2 ppb, 23.1 ppb, 104.5 ppb; 
0 3 0.45 ppb, 26.5 ppb, 128 ppb; 
PM,o, 3.8 |ig/m 3 , 43, 3183.4 ug/m 3 








CO, Ql 0.13-0.39 ppm, Q4 0.72-1.37 ppm; 


Padula et al. 






N0 2 , Ql 2.4-13.36 ppb, Q4 20.54-638.94 ppb; 


2013 [18] 






0 3 , Ql 10.49-29.05 ppb, Q4 62.65-91.92 ppb; 
PM 10 , Ql 7.9-25.24 ug/m 3 , Q4 44.09-95.32 ug/m 3 


Dadvand et al. 
2011 [29] 






Percentile 25- percentile 75: CO, 0.39-0.64 mg/m 3 ; 
N0 2 , 29.2-38.4 p.g/m 3 ; NO, 13.3-32.5 ug/m 3 ; 
0 3 , 33.2-42.4 ^g/m 3 ; PM 10 , 20.5-30.2 ^g/m 3 



Dadvand et al. 

Percentile 25-percentile 75: S0 2 , 17.6-31.2 Lig/m 

2011 [28] 





Mean: 


Mean: 




PM 10 , 28.7 ug/m 3 


PM 10 , 28.1 ug/m 3 


Marshall et al. 


N0 2 , 2.4E-2 ppm 


N0 2 , 2.4E-2 ppm 


2010 [30] 


S0 2 , 5.3E-3 ppm 


SO z , 5.1E-3ppm 




0 3 , 2.5E-2 ppm 


0 3 , 2.5E-2 ppm 




CO, 0.83 ppm 


CO, 0.85 ppm 
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Table 2. Cont. 



Studies 



Case 



Control 



Congenital Defects 



Air Pollutants 



Congenital Defects 



Air Pollutants 



Total 



Dolk etal. 
2010* [31] 



Hansen et al. 
2009 [32] 



Hwang et al. 
2008 [24] 



Gioboa et al. 
2005 [36] 



Percentile 10, median, percentile 90: 
S0 2 , 3.87 |ig/m 3 , 7.86 ug/m 3 , 14.99 ug/m 3 
N0 2 , 21.48 ug/m 3 , 35.11 p.g/m 3 , 47.78 ug/m 3 
PMio, 18.84 ug/m 3 , 21.97 ug/m 3 , 26.4 ug/m 3 
Minmum, mean, maxmum: SO2, 0, 1.5 ppb, 7.1 ppb; 
CO, 0.02 ppm, 1.1 ppm, 7.0 ppm; 
N0 2 , 1.4 ppb, 8.2 ppb, 22.7 ppb; 
0 3 , 4.3 ppb, 25.8 ppb, 54.4 ppb; 
PM10, 4.4 ug/m 3 , 18.0 ug/m 3 , 151.7 ug/m 3 



Rankin et al. 


S0 2 Q1-Q3, 




2009 [33] 


2.7-4.4 |xg/m 3 




Strickland et al. 




IQR: S0 2 , 4.0 ppb; CO, 0.3 ppm; NO z , 5.7 ppb; 


2009 * [34] 




0 3 , 29.9 ug/m 3 ; PM 10 , 14.2 ug/m 3 



Minmum, median, maxmum: 
0 3 , 16.7 ppb, 26.8 ppb, 45 ppb; 
CO, 25 pphm, 62 pphm, 277 pphm; 
NO x , 1.0 ppb, 20.2 ppb, 44.2 ppb; 
PM10, 20.8 |ig/m 3 , 57.2 ug/m 3 , 78.1 ug/m 3 
CO, Ql < 0.4 ppm, Q4 > 0.7 ppm; 
N0 2 , QK1.3 pphm, Q4 > 2.1 pphm; 
0 3 , Ql < 1.8 pphm, Q4>3.1 pphm; 
S0 2 , Ql < 1.3 ppb, Q4 > 2.7 ppb; 
PM10, Ql < 19.5 (xg/m 3 , Q4 > 29 ug/m 3 



Notes: CHD, congenital heart defects; VSD, ventricular heart defects; IQR, interquartile range.; Ql, quartile 1; Q2, quartile 2; Q3, quartile 3; Q4, quartile 4; 
* indicated cohort studies, others were designed as case-control studies; "-" no information. 
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Regarding cardiovascular anomalies, 16 combinations of pollutants -cardiac anomalies could be 
included in the meta-analysis (Figures 1-3). In all, exposure was expressed as continuous variables. 
We found a significantly increased meta-OR for exposure to NO2 and the risk of coarctation of aorta 
(OR per 10 ppb = 1.20, 95% CI (1.02, 1.41)), which is in accord with the previous meta-analysis [16]. 
For all other combinations, the combined effects were close to one and not significant. 

Figure 1. Forest plots for combinations of ventricular septal defects and pollutant 
(as a continuous variable). The size of each square represents the weight that contributes to 
the combined effect, respectively for: (A) S0 2 ; (B) PM 10 ; (C) N0 2 ; (D) CO; and (E) 0 3 . 



Study 



Glanicolo et al. 2014 
Agay-Shay et al. 2013 
Dadvand et al. 2011b 
Strickland et al. 2009 
Hansen et al. 2009 



Combined effect ^ 

(l-squared = 77.9%, p = 0.001) 



OR (95% CI) 


n 


Study 


1.230 (0.890, 1.700) 


40 


Schembari et al. 2013 


1.015(0.977, 1.055) 


493 


Agay-Shay et al. 2013 


0.964 (0.950, 0.974) 


1154 


Dadvand etal. 2011b 


0.997 (0.981,1.013) 


1654 


Hansen et al. 2009 


0.895 (0.790, 1.013) 


222 


Strickland et al. 2009 


0.986 (0.956, 1.016) 


3563 


Combined effect 



(l-squared = 58.7%, p = 0.046) 



o 



OR (95% CI) 


n 


0.653 (0.495, 0.903) 


106 


0.990 (0.940,1.040) 


493 


0.932 (0.817,1.062) 


1154 


0.706(0.455,1.103) 


222 


0.986 (0.884,1.103) 


1654 


0.928(0.842.1.023) 


3629 



.8 1 
OR per 1 pg/m 



1.2 



.8 1 1.2 
OR per 10 pg/m 3 



(A) SO2 and ventricular septal defects 

Study 



(B) PM10 and ventricular septal defects 



Schembari et al. 2013 
Agay-Shay et al. 2013 
Dadvand et al. 2011b 
Strickland et al. 2009 
Hansen et al. 2009 

Combined effect 

(l-squared = 25.9%, p = 0.249) 



<> 



OR (95% CI) 


n 


Study 


0.969 (0.793, 1.174) 


351 


Agay-Shay et al. 2013 


0.990 (0.840, 1.150) 


493 


Dadvand et al. 2011b 


0.910 (0.753, 1.098) 


1154 


Strickland et al. 2009 


1.201 (0.983, 1.418) 


1654 


Hansen et al. 2009 1 


0.768 (0.367, 1.545) 


222 


Ritz et al. 2002 


1.008 (0.923. 1.102) 


3874 


Combined effect « 



(l-squared = 91.3%, p = 0.000) 



.8 1 
OR per 10 ppb 



1.4 




OR (95% CI) 


n 


0.990 (0.650, 1.150) 


493 


3.033 (2.050, 4.485) 


1154 


0.935 (0.629, 1.330) 


1654 


0.439 (0.284, 0.661) 


222 


1.330 (1.000, 1.780) 


235 


1.109 (0.629, 1.956) 


3758 



.8 1 1.2 
OR per 1 ppm 



(C) NO2 and ventricular septal defects 

Study 
Agay-Shay et al. 2013 
Dadvand et al. 2011b 
Strickland et al. 2009 
Hansen et al. 2009 -» 
Ritz et al. 2002 

Combined effect 

(l-squared = 35.2%, p = 0.186) 



(D) CO and ventricular septal defects 



OR (95% CI) 


n 


0.950 (0.870, 1.040) 


493 


1.041 (0.951, 1.149) 


1154 


1.011 (0.967, 1.059) 


1654 


0.820 (0.660, 1.010) 


222 


1.842(0.418, 8.114) 


261 


0.999(0.963, 1.036) 


3784 



.8 1 1.4 
OR per 5 ppb 



(E) O3 and ventricular septal defects 
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Figure 2. Forest plots for combinations of atrial septal defects and pollutant (as a continuous 
variable). The size of each square represents the weight that contributes to the combined 
effect, respectively for: (A) S0 2 ; (B) PM 10 ; (C) N0 2 ; (D) CO and (E) 0 3 . 



Study 

Agay-Shay et al. 2013 
Dadvand etal. 2011b 
Strickland et al. 2009 
Hansen et al. 2009 

Combined effect 

(l-squared = 30.0%, 
p = 0.232) 



.8 



OR (95% CI) 


n 


Study 




0.988(0.948,1.026) 


534 


Schembari etal. 2013 


< 


0.984(0.957,1.013) 


274 


Agay-Shay et al. 2013 




1.000 (0.969,1.031) 


379 


Dadvand et al. 2011b 




1.182 (0.994,1.442) 


127 


Hansen et al. 2009 
Strickland et al. 2009 




0.992 (0.974, 1.011) 


1314 


Combined effect 

(l-squared = 0.0%, p = 


0.652) 



OR (95% CI) 
0.559 (0.238, 1.333) 
1.010(0.960,1.060) 



58 
534 



1.020(0.792,1.293) 274 



1.157 (0.686, 1.927) 
1.083 (0.870, 1.349) 



127 
379 



1.013(0.966,1.062) 1372 



1.4 



OR per 1 ug/m 3 

(A) SO2 and atrial septal defects 



.8 1 1.4 
OR per 10 pg/m 3 



Study 

Schembari et al. 2013 — 
Agay-Shay et al. 2013 
Dadvand et al. 2011b 
Strickland et al. 2009 
Hansen et al. 2009 

Combined effect 

(l-squared = 0.0%, p = 0.504) 



(B) PMio and atrial septal defects 



<> 



OR (95% CI) 


n 


Study 


0.909 (0.695,1.174) 


229 


Agay-Shay et al. 2013 


0.930 (0.800.1.080) 


534 


Dadvand et al. 2011b 


0.963(0.646,1.398) 


274 


Strickland et al. 2009 


1.278 (0.864,1.870) 


379 


Hansen et al. 2009 


1.240(0.537, 2.991) 


127 


Ritz et al. 2002 


0.961 (0.855,1.080) 


1543 


Combined effect 



(l-squared = 34.7%, p = 0.190) 



<> 



OR (95% CI) 


n 


0.890 (0.590,1.350) 


534 


0.868 (0.322,2.318) 


274 


0.757 (0.367,1.590) 


379 


0.733(0.451,1.208) 


127 


1.330(1.000,1.780) 


385 


0.969 (0.736,1.275) 


1699 



.8 1 1.4 
OR per 10 ppb 

(C) NO2 and atrial septal defects 

Study 
Agay-Shay et al. 2013 
Dadvand etal. 2011b 
Strickland et al. 2009 
Hansen et al. 2009 
Ritz et al. 2002 < 



.8 1 1.4 
OR per 1 ppm 



(D) CO and atrial septal defects 



Combined effect 

(l-squared = 0.0%, p = 0.706) 



OR (95% CI) 


n 


1.010(0.920,1.110) 


534 


1.062(0.868.1.280) 


274 


1.025 (0.935, 1.123) 


379 


1.030(0.770,1.370) 


127 


0.444(0.135, 1.403) 


385 


1.020(0.960.1.084) 


1699 



.8 1 1.4 
OR per 5 ppb 



(E) O3 and atrial septal defects 
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Figure 3. Forest plots for combinations of two cardiac anomalies (coarctation of aorta and 
tetralogy of fallot) and pollutant. The size of each square represents the weight that 
contributes to the combined effect. (A), (C), (E) were combined effects of coarctation of 
aorta and S0 2 , PMio and N0 2 . (B), (D) and (F) were combined effects of tetralogy of fallot 
and S0 2 , PM 10 and N0 2 . 



Study 

Dadvand et al. 2011b 
Dadvand et al. 2011a 
Dolketal. 2010 
Strickland et al. 2009 

Combined effect 

(l-squared = 67.1%. p = 0.028) 



o 



OR (95% CI) n 

0.962(0.921,1.002) 125 

0.970(0.950,1.000) 127 

1.013(0.991,1.035) 176 

1.004(0.973,1.035) 275 

0.989(0.965,1.015) 703 



.8 1 
OR per 10 ug/m 3 



1.2 



(A) S0 2 and coarctation of aorta 

Study OR (95% CI) n 



Schembarl et al. 2013 
Dadvand et al. 2011b 
Dolketal. 2010 — 



Strickland et al. 2009 

Combined effect 

(l-squared = 0.0% 



<> 



0.967(0.537,1.736) 28 

0.980(0.658,1.424) 125 

-) 1.121 (0.379, 3.289) 176 

1.103(0.884,1.380) 275 

1.062(0.886,1.272) 604 



.8 1 1.4 

OR per 10 ug/m 3 

(C) PMio and coarctation of aorta 

Study OR (95% CI) 



n 



Schembari et al. 2013 



Dadvand et al. i_ 
2011b 



Dolk etal. 2010 
Strickland et al. 2009 



Combined effect 

(l-squared = 0.0%, p = 0.563) 



1.240(1.015,1.516) 69 

0.812(0.464,1.424) 125 

1.278(0.822,1.991) 176 

1.201 (0.783,1.827) 275 

1.197(1.019,1.406) 645 



.8 1 1.4 
OR per 10 ppb 



(E) N0 2 and coarctation of aorta 



Study 

Dadvand etal. 2011b 
Dadvand etal. 2011a 
Dolket al. 2010 
Strickland et al. 2009 

Combined effect 

([-squared = 78.1%, p = 0.003) 



<> 



OR (95% CI) n 

0.945(0.906,0.985) 126 

1.000(0.980,1.020) 140 

1.029(1.006,1.054) 146 

0.985(0.954,1.015) 299 

0.993(0.964,1.023) 711 



.8 1 

OR per 10 ug/m 3 



1.2 



(B) S0 2 and tetralogy of fallot 

Study OR (95% CI) n 



Schembari et al. 2013 
Dadvand et al. 2011b 
Dolketal. 2010 
Strickland et al. 2009 

Combined effect 

(l-squared = 7.3%, 
p = 0.357) 



o 



0.704(0.350,1.459) 17 

0.651 (0.416,1.000) 126 

1.680(0.475,5.928) 146 

0.914(0.738,1.137) 299 

0.857(0.712,1.031) 588 



.8 1 1.4 
OR per 10 ug/m 3 

(D) PMio and tetralogy of fallot 

Study OR (95% CI) 



Schembari etal. 2013 
Dadvand etal. 2011b - 
Dolketal. 2010 
Strickland etal. 2009 



Combined effect 

(l-squared = 0.0%, p = 0.715) 



0.939(0.552,1.588) 49 

0.945(0.502,1.650) 126 

- 1.298(0.783,2.156) 146 

0.897(0.590,1.377) 299 

1.002(0.779.1.287) 620 



.8 1 1.4 
OR per 10 ppb 



(F) N0 2 and tetralogy of fallot 



For oro-facial clefts, seven articles were included in the meta-analysis, with one study exhibiting 
a statistically significant association [24]: the OR comparing the 4th with the 1st quartiles of the 
exposure distribution was equal to 1.40 (95% CI (1.05, 1.91)). The meta-analysis found no significant 
association for the four other pollutants (Figure 4A-D). In most cases, the studies on cleft lips, 
for which exposure was expressed on a continuous scale, were less than four, except for those 
exploring the effect of exposure to N0 2 whose result is presented in Figure 4E also showing no 
significant association. 
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Figure 4. Forest plots for a variety of pollutants and risk of orofacial defects. (A), (B), (C) 
and (D) showed combined effects of PMio, NO2, CO and O3; (E) was the combined effect 
of NO2 and cleft lip from continuous exposure risk estimates. 



Study 
Pad ula etal. 2013a 
Marshall et al. 2010 
Hwang et al. 2008 
Gllboa et al. 2005 



Combined effect 

(l-squared = 0.0%, p = 0.544) 



OR (95% CI) 
1.09(0.74,1.59) 
1.10(0.80, 1.50) 
1.00(0.82,1.21) 
-) 1.37(0.94,2.00) 



n 

75 
92 
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Study 
Padula etal. 2013a - 
Marshall et al. 2010 
Hwang etal. 2008 
Gilboa et al. 2005 



<o 



1.08(0.94,1.24) 



Combined effect 

(l-squared = 52.4%, p = 0.098) 



OR (95% CI) n 

0.69(0.47,1.02) 59 

1.30(0.90,1.80) 90 

0.90(0.58,1.37) 653 

0.82(0.56,1.19) 285 

0.92(0.76,1.14) 1087 



.8 1 1.4 

(A) PMio and cleft lip 



.8 1 1.2 

(B) N0 2 and cleft lip 



Study 

Padula et al. 2013a 1 — 

Marshall et al. 2010 
Hwang et al. 2008 
Gilboa et al. 2005 

Combined effect 

(l-squared = 72.4%, p = 0.012) 
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Hwang etal. 2008 
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Combined effect 



(l-squared = 9.7%, p = 0.345) 



.8 1 1.2 
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1.20(0.81 1.77) 
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86 
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1.09(0.70,1.69) 
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1.17(0.98,1.41) 


1117 



1.4 



(C) CO and cleft lip 

Study 



(D) 0 3 and cleft lip 



Schembari et al. 2013 
Dolk etal. 2010 
Hansen et al. 2009 
Hwang et al. 2008 



Combined effect 

(l-squared = 0.0%, p = 0.572) 



OR (95% CI) 


n 


1.126(0.894,1.410) 


130 


0.935 (0.686,1.272) 


586 


0.611 (0.647, 4.051) 


145 


0.665 (0.122,3.395) 


653 


1.066(0.892,1.275) 


1514 



.8 1 1.4 
OR per 10 ppb 

(E) N0 2 and cleft lip 



4. Discussion 



Congenital anomalies are the leading cause of neonatal mortality. The prevalence in France reported 
by the national health surveillance institute is equal to 3.3% [37]. In this systematic review, 
we assessed the effects of air pollution on the risk of congenital anomalies based on the epidemiologic 
literature. We compared 21 combinations of air pollutants and congenital anomalies types and only one 
significant result was revealed. The input data of the 21 combinations were different; it is improbable 
that the number of statistical tests performed explains the result. But, as all significant results, 
we cannot exclude that it is a chance finding. The combined risk of coarctation of aorta was found 
significantly associated with N0 2 (OR per 10 ppb = 1.20, 95% CI (1.02, 1.41)). Our meta-analysis did 
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not reveal any other significant association for cardiac anomalies and oral-facial defects, which is in 
coherence with the previous meta-analysis published in 2010 by Vrijheid et al. [16]. 

However, these authors revealed a significantly increased risk between exposure to NO2, SO2 and 
the tetralogy of fallot (OR per 10 ppb N0 2 = 1.25, 95% CI (1.02-1.51); and OR per 1 ppb S0 2 = 1.04, 
95% CI (1.00-1.08), respectively) and coarctation of aorta (OR per 10 ppb N0 2 = 1.20, 95%CI 
(1.00-1.44); OR per 1 ppb S0 2 = 1.04; 95% CI (1.00-1.08)). We did not undertake the meta-analysis 
for these combinations because the number of studies was small. In the present work, we chose to 
realize two-separate meta-analyses according to the type (continuous or categorical) of the exposure 
variable, not to introduce too much heterogeneity in the meta-estimates, whereas Vrijheid et al. 
converted all continuous to categorical variables in order to increase their statistical power. 

The present meta-analysis has several limitations. The first one deals with the sample size. 
Most of the combined effects were computed with about five studies; we did not try to compute 
combined effects for less than four studies. There was, however, a large total number of cases included 
in the meta-analysis (varying from 588 to 3874 according to the air pollutant-congenital anomaly 
combination), that provided enough statistical power to reveal significant associations. The numbers of 
cases and controls of each study are reported in Tables 3-7. We calculated the statistical power 
(a fixed alpha-risk = 5%) for the smallest sample size (Table 7) and the highest (Tables 5 and 6). 
With a statistical power equal to 90%, we will detect a statically significant health effect equal to 18% 
and 25% from the smallest and the highest sample size, respectively; note that, due to the small 
number of cases counted in the highest sample size (about 600 cases among a population of about 
1.5 million because of the cohort studies include in the meta analysis), a smaller effect could be found 
significant from the smallest sample size. The statistical power will fall to 70% and 60% 
(respectively, in the smallest and highest samples sizes) to reveal a significant increase of the risk 
equal to 10%. 

Secondly, we found differences in the methodologies adopted in the 17 individual studies included 
in the present review, which make difficult the comparisons and assessment of the combined effects. 
Certain studies relied on measurement of air quality at the birth residence [9,28-33,36], 
which is not necessarily the address during early pregnancy (the critical window time) 
and hence may lead to misclassification of exposure for several pregnant women. In addition, 
without distance-weighted calculation in the exposure assessments, the range of distance between 
maternal place and the nearest monitoring station may produce uncertain exposure estimates and thus 
lead also to exposure misclassification. 
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Table 3. Numbers of ventricular septal defects cases and number of controls in studies included in the meta analysis for 5 pollutants: 
S0 2 , PM 10 , N0 2 , CO, 0 3 . 



Numbers of cases and controls exposed under various Gianicolo et ah 
pollutants among different studies 2014 [23] 



Schembari 
etal. 2013 
[22] 



Agay-Shay 
etal. 2013 
[21] 



Dadvand 
etal. 2011 
[29] 



Strickland et al. 
2009 [34] 



Hansen et al. 
2009 [32] 



Ritz et al. 
2002 [9] 



Total Number 





case control 


case 


control 


case 


control 


case 


control 


case 


control 


case 


control 


case 


control 


case 


control 


SO, 


40 150 






493 


130,402 


1154 


4616 


1654 


713,846 


222 


1110 


235 


9049 


3798 


859,173 


PM I0 




106 


903 


493 


130,402 


1154 


4616 


1654 


713,846 


222 


1110 


235 


9049 


3864 


859,926 


NO, 




351 


2869 


493 


130,402 


1154 


4616 


1654 


713,846 


222 


1110 






3874 


852,843 


CO 








493 


130,402 


1154 


4616 


1654 


713,846 


222 


1110 






3523 


849,974 


o 3 








493 


130,402 


1154 


4616 


1654 


713,846 


222 


1110 






3523 


849,974 



Table 4. Numbers of atrial septal defects cases and number of controls in studies included in the meta analysis for 5 pollutants: SO2, PMi 0 , 
N0 2 , CO, 0 3 . 



Numbers of cases and controls exposed under various 
pollutants among different studies 


Schembari et al. 
2013 [22] 


Agay-Shay et al. 
2013 [21] 


Dadvand et al. 
2011 [29] 


Strickland et al. 
2009 [34] 


Hansen et al. 
2009 [32] 


Ritz et al. 
2002 [9] 


Total Number 




case 


control 


case 


control 


case 


control 


case 


control 


case 


control 


case control 


case 


control 


so 2 






534 


130,402 


274 


1096 


379 


715,121 


127 


635 




1314 


847,254 


PM 10 


106 


903 


534 


130,402 


274 


1096 


379 


715,121 


127 


635 




1420 


848,157 


NO, 


229 


2869 


534 


130,402 


274 


1096 


379 


715,121 


127 


635 




1543 


850,123 


CO 






534 


130,402 


274 


1096 


379 


715,121 


127 


635 


385 3000 


1699 


850,254 


0 3 






534 


130,402 


274 


1096 


379 


715,121 


127 


635 


385 3000 


1699 


850,254 
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Numbers of cases and controls exposed under various Schembari et al. Dadvand et al. Dadvand et al. Dolk et al. Strickland et al. 2009 

Total Number 

pollutants among different studies 2013 [22] 2011 [29] 2011 [28] 2010 [31] [34] 





case 


control 


case 


control 


case 


control 


case 


control 


case 


control 


case 


control 


so 2 






125 


500 


127 


508 


176 


759,817 


275 


715225 


703 


1,476,050 


PM 10 


28 


890 


125 


500 






176 


759,817 


275 


715225 


604 


1,476,432 


N0 2 


69 


2869 


125 


500 






176 


759,817 


127 


715225 


497 


1,478,411 



Table 6. Numbers of tetralogy of fallot cases and number of controls in studies included in the meta analysis for the 3 pollutants: SO2, PM10, NO2. 



Numbers of cases and controls exposed under various 
pollutants among different studies 


Schembari et al. 
2013 [22] 


Dadvand et al. 
2011 [29] 


Dadvand et al. 
2011 [28] 


Dolk et al. 
2010 [31] 


Strickland et al. 
2009 [34] 


Total Number 




case control 


case 


control 


case control 


case 


control 


case 


control 


case 


control 


S0 2 






126 


504 


140 560 


146 


759,817 


299 


715,201 


711 


1,476,082 


PM 10 


17 


890 


126 


504 




146 


759,817 


299 


715,201 


588 


1,476,412 


N0 2 


49 


2650 


126 


504 




146 


759,817 


299 


715,201 


620 


1,478,172 



Table 7. Numbers of cleft lip cases and number of controls in studies included in the meta analysis for 4 pollutants: PMi 0 , NO2, CO, O3. 



Numbers of cases and controls exposed under 

Padula et al. 2013 [18] Marshall et al. 2010 [30] Hwang et al. 2008 [24] Giloba et al. 2005 [36] Total Number 

various pollutants among different studies 





case 


control 


case 


control 


case 


control 


case 


control 


case 


control 


PM 10 


75 


200 


92 


12,925 


653 


6530 


290 


3450 


1110 


23,105 


N0 2 


59 


205 


92 


12,925 


653 


6530 


285 


3237 


1089 


22,897 


CO 


45 


157 


92 


12,925 


653 


6530 


293 


3309 


1083 


22,921 


O3 


73 


201 


92 


12,925 


653 


6530 


305 


3594 


1123 


23,250 
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Different classifications of congenital anomalies also yield some heterogeneity. Specific definition 
of defects could ease detection of significant associations. Ventricular septal defects are the most 
common congenital subgroups that were studied and a variety of classifications were found. 
For example, ventricular septal defects were classified into three subgroups (ventricular septal 
defects-permembranous, muscular and conov) in the Padula et al. study [19], according to the 
classification proposed by Botto et al. [38], whereas in others they used two subgroups (ventricular septal 
defect, muscular and permembranous) [31] or only one global group [21,29,32,36]. Another example is 
the subcategory "conotruncal defect": four studies used this outcome category to include other 
cardiovascular malformations, such as "tetralogy of fallot", "transposition of the great arteries", 
truncus arteriosus communis, "double outlet right ventricle" or "aorticopulmonary window" [9,31,32,36]; 
while, tetralogy of fallot and transposition of the great arteries were considered as an event in most 
other studies. 

Confounding factors included in the individual studies are an additional problem. Few confounders, 
such as smoking [39], parental occupation [40], maternal age [41] and season conception [42] 
have been addressed in congenital anomalies research. Season and maternal age at conception were the 
most frequent confounders considered in the studies included in the present work. Seasonal variations 
of congenital anomalies incidence have been well described, with a higher risk in summer than in 
winter [42]. The association between maternal age and the risk of non-genetic congenital anomalies is 
still unclear. Previous studies analyzing the EUROCAT database found that teenage mothers were at 
higher risk but not older mothers (35^14 years). However, among European countries, the maternal 
age patterns are less clear. This suggests the influence of social factors, exposures or living habits, 
which might be associated with maternal age [41]. The prevalence of congenital anomalies is related to 
the social status, with higher values in deprived categories [43]. Now, few studies have addressed 
socioeconomic status/deprivation as confounders or effect modifiers [22,27,28,36]. One study matched 
cases and controls according to a neighborhood socioeconomic deprivation index but without 
considering other confounders in the statistical analysis [23]. One may consider that socioeconomic 
deprivation would constitute a good approximation of well-known risk factors of congenital 
abnormalities such as smoking habits or educational level. The wide range of confounders considered 
in the individual studies included in the present review might introduce heterogeneity when combining 
the data. Because of these limitations, this meta-analysis could detect only few significant associations 
between air pollution and birth defects; however it does not mean that the hypothesis should be 
definitively disregarded. 

The different methods used in the studies, such as case and control definition, exposure assessment 
and confounding factors as all mentioned earlier, could damage the quality level of each study 
included in the meta-analysis and consequently the quality of the combined estimates. Assessing the 
quality of studies is important to understand properly each study to be used in meta analysis. It would 
be interesting to use the Newcastle-Ottawa Scale (NOS) [44] followed by quality score analysis as 
recommended by Detsky et al. [45] in order to assess the quality of each study. Then, including only 
studies with the highest quality score in the meta-analyses, we could measure more precisely the 
impact of the study quality on the point estimates. However, due to the limited number of studies in 
our meta-analyses, it was not possible to conduct this procedure. 
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Language selection may also bias the data basis. Non-English publications of relevant articles may 
have been ignored. Moreover, the risk of publication bias is inherent in systematic literature review. 
Unpublished results (probably, more likely to bear not significant results and the grey literature, 
which is not available on open sources) may distort the meta-analysis findings. So far, the two 
meta-analyses (the one published in 2010 and the present one) tend to suggest an adverse effect of air 
pollution on at least one type of birth defects, and this call for further studies in order to confirm the 
finding. We failed to assess publication bias by using funnel plots. According to the recommendations 
from Cochrane Handbook for Systematic Reviews of Interventions [46], there should be at least 
10 studies in the meta-analysis to distinguish real asymmetry. 

Air pollutants could directly exert adverse effects as pro-oxidants binding to lipid and proteins, 
therefore promoting oxidative stress and the production of free radicals, a process that may elicit a 
variety of diseases or defects [47]. This oxidative stress caused by air pollution during pregnancy has 
been pointed out in some studies [48,49]. In addition, there is recent evidence that air pollutants can 
contribute to epigenetic changes, including alteration of DNA methylation [50]. MicroRNA has been 
also studied with regard of the environmental changes and there is evidence that microRNA expression 
and regulation may be affected by environmental exposures, such as air pollution, smoking and heavy 
metal accumulation [51]. Such epigenetic modifications during pregnancy could impair normal 
embryo development and lead to birth defects. 

5. Conclusions 

Air pollution is a universal issue. Therefore, a small increase in risks may lead to serious public 
health problems. Congenital anomalies are the main causes of preterm and neonatal mortality and 
morbidity. Meta-analysis is an appropriate tool to enhance statistical power in the analysis of weak 
associations. It might shed new light on the association between air pollution and congenital anomalies 
insofar as new studies are conducted that overcome the limitations discussed in the present literature 
review. Improved exposure assessment methods, in particular more accurate spatial measurements or 
modeling, standardized definition of cases and accommodation of known or putative confounders are 
highly recommended for future congenital anomalies research on the effect of air pollution. 
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